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We report the results for nucleon decay searches via modes favored by supersymmetric grand 
unified models in Super-Kamiokande. Using f 489 days of full Super-Kamiokande-1 data, we searched 
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for p — » vK + , n — ► 1/.K" , p — > fj, + K° and p — » e + K° modes. We found no evidence for nucleon decay 
in any of these modes. We set lower limits of partial nucleon lifetime 2.3x10 , 1.3x10 , 1.3x10 
and l.OxlO 33 years at 90% confidence level for p — ► vK + , n — ► i/ff°, p — > fi + K° and p — > e + 7f° 
modes, respectively. These results give a strong constraint on supersymmetric grand unification 
models. 

PACS numbers: Valid PACS appear here 



I. INTRODUCTION 

Grand Unified Theories (GUTs) seek to unify 
the strong and electroweak forces. They are motivated 
by the apparent merging of the coupling constants of the 
strong, weak, and electromagnetic forces at a large en- 
ergy scale (~ 10 16 GeV) when low energy measurements 
are extrapolated. One of the generic predictions of GUTs 
is the instability of the proton, as well as neutrons bound 
inside the nucleus. The experimental observation of nu- 
cleon decay would provide a strong evidence of GUTs. 

In GUTs, nucleon decay can proceed via an exchange 
of a massive boson between two quarks in a proton or in 
a bound neutron. In this reaction, one quark transforms 
into a lepton and another into an anti-quark which binds 
with a spectator quark creating a meson. The favored de- 
cay mode in the prototypical GUT [2j based on an SU (5) 
symmetry ("minimal 5£/(5)") is p — > e + 7r°. For this de- 
cay, the proton lifetime scales as ~ M x , where Mx is 
the mass of the heavy vector gauge boson. In minimal 
SU(5), Mx is on the order of the coupling unification 
at 10 15 GeV/c 2 , yielding a predicted proton lifetime of 
r/B(p — > e + 7r°) ~ io 29±2 years. The first generation 
large water Cherenkov detector experiments , moti- 
vated by this prediction, observed no evidence of proton 
decay in this mode and ruled out the model. Also, the 
recent result by the Super-Kamiokande experiment ex- 
tended the previous results |5j]. It turns out that this 
contradiction of SU(5) with the experimental proton de- 
cay limit can be resolved by incorporating supersymme- 
try (SUSY) in the theories. 

Supersymmetry postulates that, for every SM particle, 
there is a corresponding "superpartner" with spin differ- 
ing by 1/2 unit from the SM particle |6j. The additional 
particles stabilize the renormalization of the Higgs boson 
and address the so-called "hierarchy problem" . When 
one incorporates the superpartners into the calculation 
of the running of the coupling constants, the convergence 
of the coupling constants occurs at an unification scale 
about one order of magnitude larger than that of minimal 
SU(5). Since the proton decay rate via p — > e + ir° scales 
as M^ 4 , this leads to a suppression of about four orders 
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of magnitude in the rate, consistent with experimental 
non-observation of p — > e + 7r°. Furthermore, while in the 
minimal SU(5) model, the three coupling constants do 
not quite meet at a single point within three standard 
deviations, they meet together at a single point in the 
minimal SUSY SU(5) model 0. 

However, in many SUSY GUT models, other dominant 
nucleon decay modes occur via dimension five operator 
interactions with the exchange of a heavy supersymmet- 
ric color triplet Higgsino These interactions suppress 
transitions from one quark family in the initial state to 
the same family in the final state. Since the only second 
or third generation quark which is kinematically allowed 
is the strange quark, an anti-strange quark typically ap- 
pears in the final state for these interactions. The anti- 
strange quark binds with a spectator quark to form a K 
meson in the final state. Thus, SUSY GUTs favor nu- 
cleon decays in p — > vK + and n — > DK° modes. The 
predictions for the nucleon lifetime in SUSY GUT mod- 
els, however, varies widely, and may even be suppressed, 
since there are many new free parameters introduced be- 
cause of the supersymmetry breaking. 

In the minimal SUSY 5*17(5) GUTs, the partial proton 
lifetime is estimated to be r/B(p -> vK + ) < 2.9 x 10 30 
years 9] . Other models have been proposed beyond those 
based on SU(5). In particular, models based on SO(10) 
symmetry have become popular in light evidence for neu- 
trino mass |10|, which they naturally accommodate. An 
important property of the 5*0(10) symmetry is that there 
is a heavy right-handed neutrino in a multiplet contain- 
ing the matter fields. In addition, all matter fields of 
one generation can be contained in a single multiplet, in 
contrast to theories based on SU (5) where they must be 
broken into two separate representations. 

One class of models uj. predicts the partial pro- 
ton lifetime for p — > vK + mode to be less than 10 34 
years, which is within the observable range of Super- 
Kamiokande. In addition, the same mechanism which 
gives mass to the neutrinos provides a new set of di- 
mension five operators through which the proton can de- 
cay. A consequence of this is that the prediction for 
p — » fi + K° decay rate is within a factor of 10 of the 
p — > vK + decay rate. 

These decay modes, favored by SUSY GUT models, 
have been searched for in water Cherenkov detector 
0, and iron calorimeter [l2L fl3| experiments. The 
best limits on the partial nucleon lifetime via p — + vK + , 
n -» i?K , p -> n + K°, and p -v e + K° are 6.7 x 10 32 
years 8.6 x 10 31 years 0, 1.2 x 10 31 years 0HE!, 
and 1.5 x 10 31 years 0, respectively. 
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II. SUPER-KAMIOKANDE DETECTOR 

The Super-Kamiokande detector ^4| is a 50 kton wa- 
ter Cherenkov detector located at the Kamioka Ob- 
servatory of Institute for Cosmic Ray Research in the 
Kamioka mine. It lies 1,000 m underneath the top of Mt. 
Ikenoyama, (i.e. 2,700 m water equivalent underground), 
resulting in a cosmic ray muon rate of 2.2 Hz, a reduc- 
tion of 10~ 5 compared to the rate at the surface. The 
detector is optically separated into two regions, the inner 
and outer detectors (ID and OD). The ID of the Super- 
Kamiokande-I detector, which operated from April 1996 
to July 2001, was instrumented with 11,146 50-cm diam- 
eter inward facing photomultiplier tubes (PMTs) which 
provide 40% photocathode coverage. This photocathode 
coverage makes it possible to detect low energy electrons 
down to ~5 MeV. The ID volume is the sensitive re- 
gion for nucleon decay searches and has a total fiducial 
volume of 22.5 kton, defined as a cylindrical volume with 
surfaces 2 m away from the ID PMT plane. The OD com- 
pletely surrounds the ID and is instrumented with 1,885 
20-cm diameter outward facing PMTs equipped with 60 
cm x 60 cm wavelength shifter plates to increase light 
collection efficiency. The main purpose of the OD is to 
tag incoming cosmic ray muons and exiting muons in- 
duced by atmospheric neutrinos. A detailed description 
of the S up er- K amiokande-I detector can be found else- 
where [l4j . 



III. DATA SET 

Data used for this analysis is taken during the period 
from May 1996 to July 2001. It corresponds to 1489 days 
of data taking and an exposure of 92 kt-year. 



IV. NUCLEON DECAY EVENT AND 
BACKGROUND SIMULATION 

In order to search for nucleon decay in the Super- 
Kamiokande detector, it is critical to understand the 
signal event signature and background characteristics. 
We simulate specific decay modes of nucleon decay as 
well as background events using a signal event gener- 
ator and the standard Super-Kamiokande atmospheric 
Monte Carlo (MC). By comparing the characteristics of 
these signal and background simulated events in detail, 
we establish the optimum event selection criteria. When 
we limit the nucleon decay events to occur only in the 
Super-Kamiokande fiducial volume, the only significant 
background to nucleon decays originate from atmospheric 
neutrino interactions. Once the selection criteria are es- 
tablished, the detection efficiency is then estimated by 
analyzing the nucleon decay MC sample, and the ex- 
pected background is estimated by analyzing the atmo- 
spheric neutrino MC sample. 



A. Nucleon decay event simulation 

Nucleon decay in water can occur from a free proton 
in the hydrogen nucleus or from a bound nucleon in the 
oxygen nucleus. It is relatively simple to simulate the free 
proton decay using the Super-Kamiokande MC. However, 
simulation of the bound nucleon decays requires special 
care, because of various nuclear effects experienced by 
the daughter particles from nucleon decay before they 
exit the nucleus. 

Nucleons bound in oxygen have Fermi momentum and 
nuclear binding energy. In our simulation, we use the 
Fermi momentum distribution measured by an electron- 
ic scattering experiment |l5j. For nucleon decays in an 
oxygen nucleus, the nucleon mass must be modified by 
nuclear binding energy. The modified nucleon mass m' N 
is calculated by m! N = - E^nd where is the nu- 
cleon rest mass and Eund is the nuclear binding energy. 
Yamazaki and Akaishi 0] estimated the effective nu- 
cleon mass when the nucleon decays in 16 O. Ten percent 
of decays are from a nucleus which has wave functions 
correlated with other nucleons in the nucleus. Figure ^ 
shows the invariant proton mass distribution in 16 O used 
in the p — * vK + MC simulation. The correlated de- 
cays produce the broad nucleon mass distribution below 
around 850 MeV/c 2 . 




1 r 
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invariant proton mass in 16 0(MeV/c 2 ) 

FIG. 1: The invariant proton mass distribution in 16 in the 
p — » vK + MC simulation jig. The distributions of s- and 
p-states are expressed by Gaussian function G(mean, RMS). 
It is G(938.3 - 39.0, 10.2)) MeV/c 2 for s-state and G(938.3 - 
15.5,3.82)) MeV/c 2 for p-state. 

When a nucleon decays in the oxygen nucleus, the re- 
maining nucleus can be in an excited state and it can emit 
prompt gamma-rays during its subsequent de-excitation 
process. This process was studied by Ejiri |17| and we 
use his results in this analysis. This prompt gamma-ray 
provides us with a powerful tool to tag p — * PK + events. 
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The decay products of a nucleon bound in oxygen can 
interact hadronically with the remaining protons and 
neutrons in the residual 15 N nucleus before exiting. The 
initial position of the nucleon is assumed to have the 
Woods-Saxon density distribution [l8| : 



p{r) = 



P(0) 



l + exp(2 



r — a ^ 
b t 



(i) 



where a (=1.07 A 1/3 = 2.69 fin for 16 0) is the average 
of nuclear radius, 2b (=0.82 fm) is the thickness of the 
nuclear surface, and r is the distance from the center 
of the nucleus. Interactions relevant to this analysis are 
K+N and K°N elastic scatterings, and K+N and K°N 
inelastic scatterings via charge exchange. The K+p re- 
actions can proceed elastically or inelastically. Inelastic 
interactions are due to the K+p — > K+A. A partial wave 
analysis for the scattering amplitudes was performed by 
Hyslop et al. 01 by doing a global fit to many data sam- 
ples. The maximum momentum for the kaon from the 
nucleon decays is about 600 MeV/c. Below 800 MeV/c, 
K+p scattering has an extremely small contribution from 
inelastic scattering. Therefore a K+ from a proton decay 
via p — > vK+ effectively experiences only elastic scatter- 
ings with protons in the residual nucleus. Since p — > vK+ 
are identified through the detection of the daughter par- 
ticles of the kaon decay at rest, K + p scattering does not 
affect the detection efficiency of p — > vK+ . The charge 
exchange reaction K+n — > K°p can reduce the efficiency 
for detecting p — * vK+ events. It is important to esti- 
mate what fraction of K+ is lost due to this effect. The 
reaction was measured for low momentum kaons (250 
to 600 MeV/c) by Glasser et al [2(j yielding cross sec- 
tions ranging from 2.0 ± 0.18 mb at K+ momentum of 
250 MeV/c to 6.4 ±0.56 mb at K+ momentum of 590 
MeV/c. To estimate the fraction of K+ lost due to this 
effect, a MC simulation is performed. K+ are started at 
random points in the nucleus according to the Woods- 
Saxon density distribution (Equation^). If there is an 
interaction, Pauli blocking is taken into account by re- 
quiring the momentum of the recoil nucleon to be above 
the Fermi surface momentum {pf)- 
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(2) 



where p(r) is the same as defined in Equation ^ From 
this simulation, it is estimated that 1% of K+ from 
p — » vK+ decays in total are lost because of this charge 
exchange reaction. From isospin symmetry, the K°N re- 
actions have essentially the same magnitude as the K+N 
reactions. 

We simulate propagation of the produced particles 
and Cherenkov light in water by custom code based on 
GEANT [23. The propagation of charged pions in wa- 
ter is simulated by custom code based on for less 
than 500 MeV/c and by CALOR for more than 500 
MeV/c. 



In the n — > DK°, p — > p+K° and p — > e+K° searches, 
only decays to K s are studied because the lifetime of the 
K® is so long and many of them interact in water before 
decaying. The effect of K° regeneration is small in K s 
decay searches. 

In total, 50,000 MC events are generated for each de- 
cay mode to find the signature of nucleon decays and to 
estimate detection efficiencies. Of these, 34,664, 34,561, 
34,648 and 34,572 events are in the fiducial volume in the 
p -> uK+, n -> DK°, p -> p+K° and p -> e+K° MCs, 
respectively. 



B. Atmospheric neutrino background 

The most significant background to nucleon decay 
comes from the atmospheric neutrino interactions in the 
detector. Atmospheric neutrinos are produced in colli- 
sions of cosmic rays with air molecules in the atmosphere 
of the earth. Primary cosmic rays, mostly protons, inter- 
act hadronically with air molecules creating tt+ , K+, K°, 
and other mesons. Neutrinos are then produced from the 
chain decay of these mesons. The production of atmo- 
spheric neutrinos has been calculated in great detail by 
many authors. For this analysis, we use the calculation 
by Honda et al [25|. 

Because the cross section for neutrinos to interact with 
matter is extremely small, they can travel unscattered 
through the earth and interact with a nucleon in the wa- 
ter of Supcr-Kamiokande via the weak interaction. A 
generic interaction is 



+ N -> I + N' + X, 



(3) 



where N and N' are the initial and final state nucleons, 
I is the outgoing lepton associated with v\ , and X can be 
other possible hadronic particles such as pions. Because 
some of these interactions result in topologies similar to 
those of nucleon decays, they present a challenging back- 
ground to nucleon decay searches. 

We use the same atmospheric neutrino MC simulation 
|20| that is used for Super-Kamiokande neutrino oscil- 
lation studies with slight modifications. The modifica- 
tions made for this analysis are as following: for single 
pion production, the so-called axial vector mass which 
appears in the neutrino-nucleon differential cross section 
is set to be 1.2 GeV/c 2 to estimate conservatively. For 
deep inelastic scattering, we use a model motivated by 
Bodek and Yang |27| . 

In order to estimate the atmospheric neutrino back- 
ground events for the nucleon decays involving kaons, 
the production of kaons through the baryon resonances 
is included based on the model of Rein and Sehgal [2^| . 
The cross section for kaon production is one or more than 
one order of magnitude smaller than that of single pion 
production. In addition, many charged kaons produced 
in neutrino interactions emit Cherenkov light while kaons 
from proton decay do not, and K production is accompa- 
nied by a A baryon which decays into either pir~ or nir° 
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in most cases. Therefore these background events can be 
distinguished from nucleon decay events, and their con- 
tributions to the background can be safely ignored. 

In order to estimate the backgrounds from the atmo- 
spheric neutrinos to nucleon decays at a high precision, 
we use a total of 100 year equivalent sample of atmo- 
spheric neutrino MC simulation events for this analysis. 



V. EVENT SELECTION AND 
RECONSTRUCTION 

We apply the first stage reduction for the data to re- 
move major background. Then we reconstruct physical 
quantities for data and also signal and background MC 
events. Data reduction and event reconstruction is the 
same as for the neutrino oscillation analysis |2fij . 

A. Event Selection 

The trigger threshold for events used in this analysis is 
set to 29 PMT hits corresponding to an electron equiva- 
lent energy of 5.7 MeV. The trigger rate at this threshold 
is about 10 Hz. The majority of the collected data are 
the events originating from cosmic muons and radioac- 
tivities. An event selection process is applied to remove 
these events. Almost all cosmic muons are rejected re- 
quiring no significant OD activities. After applying event 
reconstruction, events are rejected when their electron 
equivalence energy is less than 30 MeV corresponding 
to 197 MeV/c for muons or when the event vertex is 
located outside of the fiducial volume. The remaining 
event rate after this first stage reduction is about eight 
events/day. Almost all of these events originate from at- 
mospheric neutrino interaction. 



or pion particle assumption. The momentum of a pion 
is determined by Cherenkov angle as well as the sum of 
PEs. The momentum resolution for a 236 MeV/c muon 
from p — > vK + ,K + — » is 3% and mass resolution 

for tt° from p -> DK+,K+ -> tt+tt is 9%. Figure |2 
shows the reconstructed tt° invariant mass distribution 
in the p -> DK+,K+ -> tt+tt MC. The detection effi- 
ciency of Michel electrons is estimated to be 80% for //+ 
and 63% for fjr. 




invariant mass (MeV/c 2 ) 

FIG. 2: The reconstructed 7r° invariant mass distribution 
in the p — > vK + ,K + — ► tt+7t° MC after applying selection 
criteria (B1,B2) (See Section IvTA) . 



B. Event Reconstruction 

Physical quantities are reconstructed for the events re- 
maining from the event selection process. In the recon- 
struction, certain quantities are determined such as the 
event vertex, the number of visible Cherenkov rings, the 
direction of each ring, particle identification, momentum 
and number of Michel electrons. The vertex resolution is 
34 cm for electrons and 25 cm for muons in single ring 
events. In the p — > DK + , K + — * /ti+i^ and K + — > tt + tt° 
MC, the vertex resolution is 47 and 37 cm, respectively. 
In the n -» vK°, p -> ^i+K" and p -> e+K° MC, bet- 
ter resolutions are obtained. Each ring is identified as 
e-like (e ,7) or //-like (/x+,tt+, proton) based on a likeli- 
hood analysis of Cherenkov ring pattern and Cherenkov 
angle. In a single ring electron or muon event, the mis- 
identification probability is 0.5%. Momentum is deter- 
mined by the sum of photo electrons (PEs) after correct- 
ing for light attenuation in water, PMT angular accep- 
tance, PMT coverage, and the assigned electron, muon, 



C. Calibration 

An important aspect of nucleon decay analyses is the 
choice of selection criteria. These criteria are dependent 
upon the energies of the particles produced in the vari- 
ous interactions. For this reason, the energy deposition 
and detector response of particles traversing the detec- 
tor must be accurately modeled in the MC simulation. 
Energy scale calibrations are performed by comparing 
the data and MC sample of Michel electrons from stop- 
ping cosmic ray muons, 7r°s produced in atmospheric neu- 
trino interactions and stopping cosmic ray muons. The 
Cherenkov angle and the range are used to obtain muon 
momentum for stopping muons. These samples provide 
energy calibration varying from several tens MeV (Michel 
electrons) to several GeV (stopping cosmic ray muons), 
entirely covering the energy range relevant to nucleon de- 
cay searches. Based on these samples, the energy scale is 
reproduced in the MC within 2.5%. 
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VI. ANALYSIS 

During 1,489 days of livetime corresponding to 92 
kton-year of exposure, 12,179 events are recorded after 
the first stage reduction process. Of these, 8,207 are iden- 
tified as single ring. 

It has been well known for some time that there is a 
deficit of events induced by atmospheric muon neutrino 
oscillations I n order to estimate accurately the ac- 
tual number of background events, the atmospheric neu- 
trino flux is normalized based on the observed deficit of 
/j,-like events due to well established neutrino oscillation 
phenomena and the observed number of e-like events. 
The normalizing factor is determined as multiplying 0.67 
for charged current events and 1.07 for charged current 
v e and all neutral current events. Finally we normalize 
the 100 year MC sample (2,246 kton-years) to the 92 
kton-year sample of Super-Kamiokande exposure. Based 
on these simulated signal and background events, we de- 
termine selection criteria for each nucleon decay mode, 
which is discussed in the following sections. All detection 
efficiencies given below include all relevant branching ra- 
tios. 

A. Search for p — > i>K + 

This is the primary mode of proton decay in SUSY 
GUT models. A search for p — > vK + was described in 
our previous paper [^j based on the 45 kton-year data. 
In that paper, we described three methods to search for 
this mode using the two dominant decay modes of the 
kaon, K+ -> tt+tt (Method 1) and K+ -> fJ, + v^. For 
K + — > /i + ^ t , we looked for both the prompt gamma- 
ray from the de-excitation of the residual excited 15 N 
nucleus (Method 2: prompt gamma-ray search) and for 
an excess in the momentum distribution of /Lt-like events 
at p^ — 236 MeV/c (Method 3: mono-energetic muon 
search). In this paper, we present proton lifetime limits 
derived using these methods including improved analyses 
which increase the detection efficiency and reduce the 
background. 

1. Method 1: K + -> 7r + 7r° search 

If a proton decays to v K + , the K + has low enough mo- 
mentum that the majority stop before decaying. There- 
fore, when the K + decays to a 7r + and a 7r°, these two 
particles go back-to-back and the 7r° momentum is ex- 
pected to be mono-energetic at 205 MeV/c. To de- 
tect this type of events, the following criteria are re- 
quired: (Al) two e-like rings, (A2) one Michel elec- 
tron, (A3) 175 MeV/c < p n o < 250 MeV/c, (A4) 85 
MeV/c 2 < m w o < 185 MeV/c 2 , (A5) 40 PE < Q w + < 100 
PE, (A6) Qres < 70 PE. The reconstructed total momen- 
tum and invariant mass consistent with the ir° from the 
two e-like rings are defined as p^o and m^o , respectively. 



The 7r + momentum is so close to the Cherenkov thresh- 
old that the Cherenkov ring is not detected in most cases. 
However, since many backgrounds survive criteria (Al- 
A4), we use the Cherenkov light produced by the 7r + in 
addition. Q„+ is sum of PEs corrected or light atten- 
uation and PMT acceptance, which is observed in the 
PMTs within 40° half opening angle opposite to the ir° 
direction. Q res is sum of PEs in the remaining PMTs 
after rejecting the area within the 90° half opening an- 
gle toward the two gamma-ray directions and the Q„+ 
searched area. We use the criterion (A6) for background 
rejection. The total detection efficiency for this method 
is estimated to be 6.0% and the background is estimated 
to be 0.6 events. Figure shows the distribution of the 
Qn+ versus p w o for data, the atmospheric neutrino MC 
and p — •> DK + MC. The main sources of background are 
events by single pion production and deep inelastic scat- 
tering. The incoming neutrino energy of the background 
is typically between 0.6 GeV and 2 GeV. In the data, no 
events pass the selection criteria (A1-A6). 



2. Method 2: K + —* fi + Vfi, prompt gamma-ray search 

When a proton decays in one of the inner shells of the 
16 nucleus, the residual 15 N nucleus is left in an ex- 
cited state. This state quickly de-excites with a certain 
probability of emitting gamma-rays. The most proba- 
ble residual state is a p 3 / 2 6.3 MeV state which leads 
to an emission of single 6.3 MeV gamma-ray [T^ . Since 
the K + is below the Cherenkov threshold and the K + 
lifetime is 12.4 ns, we can separate the de-excitation 
gamma-ray from the /x + signal. The /i + momentum is 
mono-energetic (236 MeV/c), because the K+ decays at 
rest. By requiring a prompt gamma-ray signal as well as 
the mono-energetic muon and an electron from the muon 
decay, most backgrounds are eliminated. 

In order to search for the prompt gamma-ray, three 
quantities must be defined. The first is t^ which is a ref- 
erence time associated with the detection of the muon. 
The second is to which is the time to begin a backward 
search in time of flight (TOF) subtracted timing distri- 
bution for the earlier hits from the prompt gamma-ray. 
Finally, t 7 is the time associated with the detection of the 
gamma-ray. The reference point t^ corresponding to the 
muon is found by searching for the point in time when 
AATfrjt/At is maximum where Nhu is number of PMT 
hits. The starting time to is defined as the first point 
less than i M where dNhit/dt = 0. In the t 1 search, PMTs 
which are within a cone with a 50° half opening angle 
with respect to the muon are removed. Removing these 
tubes enables the search for the prompt gamma-ray to 
start at a closer time to t^. A 12 ns timing window is 
slid backward starting with its trailing edge at to. The 
values t 7 and Nhu-y are determined by maximizing the 
number of hits in the 12 ns sliding window. N^itj is 
the maximum number of hits at t 7 in the 12 ns sliding 
window. 



7 





140 






120 






100 












80 




+ 






o 


60 






40 






20 






atm v MC 



p^vKt MC 




200 



400 



200 
p^o (MeV/c) 



400 



200 



400 



FIG. 3: The distributions of the versus p n o for events that satisfy criteria (A1,A2,A4,A6). Each figure shows: (left) data, 
(middle) atmospheric neutrino MC and (right) p — > vK + MC. The box shows criteria (A3,A5). 



Using these quantities, the following selection criteria 
are applied: (Bl) one /z-like ring, (B2) one Michel elec- 
tron, (B3) 215 MeV/c < p p < 260 MeV/c, (B4) pro- 
ton rejection, (B5) t^ — t^ < 100 ns, (B6) 7 < Nhitj < 
60. Criterion (B4) is applied for rejecting backgrounds 
caused by poor vertex reconstruction. Most of this type 
of backgrounds is recoil protons produced by neutral cur- 
rent interactions. Since the particle type is assumed as 
a muon in the vertex fit, the reconstructed vertex posi- 
tion of a proton is not accurate. Therefore a fake peak, 
which mimics gamma-rays, is sometimes produced in the 
TOF subtracted timing distribution. In order to reject 
these events, two cuts, g > 0.6 and d^ e < 200 cm, are 
applied, where g is a goodness of the TOF subtracted 
timing distribution and d^ ie is the distance between the 
Michel positron vertex and the muon stopping point. In 
proton decay events, d^e should be close to zero. 

Passing the p — > DK + and 100 year atmospheric neu- 
trino MC events through criteria (B1-B6), the detec- 
tion efficiency and background are estimated to be 8.6% 
and 0.7 events, respectively. The efficiency includes all 
branching ratios. Figure 0] shows N^it-y distribution for 
the 100 year sample of atmospheric neutrino MC normal- 
ized by livetime and neutrino oscillation, p — > PK + MC, 
and 92 kton-yr sample of data. No events survive these 
selection criteria. 



3. Method 3: K + — » ^t + *-V, mono- energetic muon search 

We search for a 236 MeV/c mono-energetic muon us- 
ing events which are not selected in Method 2. To tag 
the muon, criteria (B1,B2) are applied, which are defined 
as in Method 2. Moreover, to obtain an independent 
event sample from Method 2, we apply criterion (B7) 
Nhit-y < 7. Figure shows the muon momentum dis- 
tribution of data compared with best fitted atmospheric 
neutrino MC. No significant excess is observed in the sig- 




FIG. 4: The Nhitj distributions for events that satisfy crite- 
ria (B1-B5). The full circles show data and the dashed line 
shows atmospheric neutrino MC normalized by livetime and 
neutrino oscillation. The solid line shows p — > vK + MC. The 
arrows show criterion (B6), where proton decay candidates 
with a prompt gamma-ray are expected. 



nal region (B3). In this region, there are 181 events with 
a best fitted background of 200 events. 



B. Search for n — > vK 

Many SUSY models also predict the nucleon decay 
mode, n — > DK°. We search for n — » DK° using the 
Kg — » 7r°7r° and K s — > tt + it~ decay chain. 

For the K s — > 7r°7r° search, the following criteria are 
applied: (CI) three or four e-like rings, (C2) zero Michel 
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muon momentum distribution for events that satisfy criteria 
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electrons, (C3) 200 MeV/c < p K o < 500 MeV/c, (C4) 
400 MeV/c 2 < m K o < 600 MeV/c 2 , where p K <> and m K o 
are total momentum and invariant mass, assuming the 
decay sequence of K°. Because it is difficult to identify 
four showering Cherenkov rings, many background events 
remain. About half of the background events come from 
deep inelastic scattering. Figure El shows px° versus m^o 
distributions after applying criteria (C1,C2). The criteria 
select 14 events in the data with an expected background 
of 19 and a selection efficiency of 6.9%. 

For the Kg — > ir + Tr~ search, the following criteria are 
applied: (Dl) two //-like rings, (D2) zero or one Michel 
electron, (D3) 200 MeV/c < p K o < 500 MeV/c, (D4) 450 
MeV/c 2 < m K o < 550 MeV/c 2 . Since a large fraction 



of 7r or 7r~ momenta are below the Cherenkov thresh- 
old, the efficiency of finding both rings is low. Figure 
shows pk o versus ttlk ° distributions after applying crite- 
ria (D1,D2). Twenty events are observed in the data with 
the detection efficiency of 5.5%, while 11.2 background 
events are expected. Most of the background events are 
produced by charged current single pion production. 



C. Search for p 



In a SUSY 5O(10) model HH with neutrino mass, 
p — > ji + K is an important decay mode. We search for 
p — > n + K° using Kg — > 7r°7r° and K s — * tt + it~ decay 
chain. In this mode, the total invariant proton mass and 
momentum can be reconstructed and the backgrounds 
can be significantly reduced. 



applied: (El) 2-4 e-like rings and one //-like ring, (E2) 
zero or one Michel electron, (E3) 400 MeV/c 2 < m K o < 
600 MeV/c 2 , (E4) 150 MeV/c < Pf , < 400 MeV/c, 
(E5) p p < 300 MeV/c, (E6) 750 MeV/c 2 < m p < 1000 
MeV/c 2 , where m p and p p are the invariant mass and 
total momentum assuming the decay sequence of proton, 
respectively. The momentum of the // + , p^, is deter- 
mined using the //-like ring. Figure [S] shows p p versus m p 
distributions after applying criteria (E1-E4). No events 
survive in the data with the estimated 5.4% detection 
efficiency and 0.4 events of the expected background. 

For the K s — » search, two different methods are 

applied, because ring finding efficiency of the two pions 
are low. The first method (Method 1) uses the following 
criteria: (Fl) two //-like rings, (F2) two Michel electrons, 
(F3) 250 MeV/c < P p < 400 MeV/c, (F4) p p < 300 
MeV/c. We assume the more energetic ring as muon 
and the other as charged pion. Figure [5] shows p p versus 
p^ distributions after applying criteria (F1,F2). The de- 
tection efficiency and the expected background are 7.0% 
and 3.2 events, respectively; 3 events are observed in the 
data. The second method (Method 2) requires: (Gl) 
three rings, (G2) one or two Michel electrons, (G3) 450 
MeV/c 2 < m K o < 550 MeV/c 2 , (G4) p p < 300 MeV/c, 
(G5) 750 MeV/c 2 < m p < 1000 MeV/c 2 . Figure EH 
shows the p p versus m p distributions after applying cri- 
teria (G1-G3). The detection efficiency and the expected 
background are estimated to be 2.8% and 0.3 events, re- 
spectively Most of the background events come from 
charged current v p single pion production. No surviving 
events are observed in the data. 



D. Search for p -» e + K° 

One of the supersymmetric theories, based on flavor 
group (S3) 3 predicts that p — > e + K° occurs at a compa- 
rable rate with p — > vK + and n — > PK° [33. For the 
search for p — > <= + 



e 1 K° mode, three methods are applied in 



the same way as p — * [i + K° search. 



For the K s 



For the K° s 



search, the following criteria are 



search, the selection criteria are: (HI) 3-5 e-like rings, 
(H2) zero Michel electrons, (H3) p p < 300 MeV/c, (H4) 
750 MeV/c 2 < m p < 1000 MeV/c 2 . We assume the most 
energetic ring as electron and others as gamma-rays from 
the 7T° decays. Figure ^2 shows p p versus m p distribu- 
tions after applying criteria (H1,H2). The detection ef- 
ficiency and the expected background are 9.2% and 1.1 
events, respectively. One candidate event is observed in 
the data. 

For the K s — > tt + tt~ search, two different methods are 
applied. The first method (Method 1) uses the following 
criteria: (II) one //-like ring and one e-like ring, (12) one 
Michel electron, (13) 250 MeV/c < p e < 400 MeV/c, (14) 
p p < 300 MeV/c, where p e is the electron momentum de- 
termined using the e-like ring. Figure IT21 shows p p versus 
p e distributions after applying criteria (11,12). The de- 
tection efficiency and the expected background are 7.9% 
and 3.6 events, respectively. Most of the background 
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events are produced by charged current v e single pion 
events. Five events are observed in the data. The second 
method (Method 2) uses the criteria: (Jl) three rings 
(at least one e-like), (J2) zero or one Michel electron, 
(J3) 450 MeV/c 2 < m K o < 550 MeV/c 2 , (J4) p p < 300 
MeV/c, (J5) 750 MeV/c 2 < m p < 1000 MeV/c 2 . Fig- 
ure El shows p p versus m p distributions after applying 
criteria (J1-J3). The detection efficiency is estimated to 
be 1.3% with 0.04 expected backgrounds. No surviving 
events are observed in the data. 



E. Systematic uncertainties 

In the detection efficiency, we consider the following 
common systematic uncertainties in every search: imper- 
fect knowledge of light scattering in water, energy scale, 
and particle identification. The uncertainty of the rate 



of light scattering in water is estimated to be 20%, which 
mainly effects on the ring finding efficiency. The energy 
scale uncertainty is estimated to be 2.5% by the calibra- 
tion described in Section lY CI For the modes which have 
a charged pion from kaon decay, we consider the imper- 
fect knowledge of a charged pion-nucleon cross section in 
water which is estimated to be 10% by comparing with 
our MC and experimental data j^. For the p — * n + K° 
and p — > e + K° search, we consider the uncertainty of 
the Fermi momentum which is estimated to be about 5% 
from model differences [l5j| . In addition, we consider the 
uncertainty of de-excitation gamma-ray emission proba- 
bilities in the p — > DK + , prompt gamma-ray search. It 
is estimated to be 15% for 6.3 MeV gamma-ray and 30% 
for other gamma-rays |17| . From these sources of sys- 
tematic uncertainty, the total contribution to the uncer- 
tainty of the detection efficiency in the p — > i>K + , prompt 
gamma-ray search is estimated to be 20%. The main 
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FIG. 11: The distributions of the total momentum versus invariant mass for events that satisfy criteria (H1,H2). Each figure 
shows: (left) data, (middle) atmospheric neutrino MC and (right) p — > e + K° MC. The box shows criteria (H3,H4). 
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uncertainty comes from the uncertainty of de-excitation 
gamma-ray emission probabilities. The total contribu- 
tions in the p — > DK + , mono-energetic and K + — > 7r + 7r 
search are estimated to be 2.5% and 8.8%, respectively. 
Most of the uncertainties in the K + — » 7r + 7r° search come 
from the imperfect knowledge of charged pion-nucleon 
cross section in water and light scattering in water, which 
are estimated to be 4.8% and 6.7%, respectively. In the 
n — > DK° search, the total contributions are estimated 
to be 16% and 14% for the K% 



tt°tt° and K9 



search, respectively. The main source of uncertainty 
comes from the uncertainty in light scattering in water. 
In all methods of the p — * /i + K° and p — > e + K° searches, 
the total contributions are less than 20%. 

In the background estimation, we consider the follow- 
ing common systematic uncertainties in all searches: im- 
perfect knowledge of atmospheric neutrino flux, neutrino 
cross sections, energy scale, and particle identification. 
Because we use the scaled MC, the absolute normaliza- 
tion error becomes smaller. However, the background 
uncertainty from the atmospheric neutrino flux normal- 
ization is conservatively estimated to be 20% j^B 1 ]. The 



imperfect knowledge of the cross sections are considered 
to be 30% for quasi-elastic and elastic scattering, 30% for 
single meson production and 50% for deep inelastic scat- 
tering. From these sources of systematic uncertainty, the 
total contributions to the uncertainty of the background 



in the n -> DK°, K% -> 7r°7r° and K° s 



search arc 



estimated to be 44% and 41%, respectively. The contri- 
butions in many other modes become more than 50%, be- 
cause background statistics of the atmospheric neutrino 
MC are poor. The uncertainties in all of the searches are 
summarized at Table [i] 



F. Results 

In the absence of any significant nuleon decay signa- 
ture, we interpret our results as lower limits of nucleon 
partial lifetime for each decay mode using the following 
method |31j . 

First, based on Bayes theorem, we calculate the nu- 
cleon decay probability, P(T\m), as follows: 



J 



P{T\n t 



A 



-(rAiei+fci) 



(rAiCi 



r 



-P(r)P(A I )P(e,)P(6 I )dA l cie I d6 i . 



(4) 



Here a Poisson distribution is assumed for the nucleon 
decay probability; m is the number of candidate events 
in the i-th nucleon decay search; T is the total decay 
rate; A* is the corresponding detector exposure; is the 
detection efficiency including the meson branching ratio; 
and bi is the expected background. In our search, Aj is 
3.05xl0 34 proton-year for proton decay and 2.44xl0 34 
neutron-year for neutron decay. P(r) is the decay rate 
probability density. We assume P(T) as one for T > 
and otherwise zero. The uncertainties of detector expo- 
sure (P(Aj)), detection efficiency (P(e<)) and background 
(P(bi)) are expressed as follows: 



P(Ai) = S(Xi - Ao,i) 



(5) 



P(e i ) = e -( £ *- e o.*) 2 /2<i (0 < e 4 < 1, otherwise 0) (6) 



P(bi) = - 

b. 



00 e" 6 ' (6')""'* 



db' 



i JO 



rib,, 



(0 < h, otherwise 0), (7) 



where Ao.i is the estimated exposure, eo,j is the estimated 
detection efficiency, a e ^i is the estimated uncertainty in 
the detection efficiency, n^j is the number of background 
events, C, is the MC oversampling factor, and (Jh,\ is the 



uncertainty in the background. Because the uncertainty 
in the exposure is small, a delta function is assumed for 
P(Ai). The lower limit of the nucleon decay rate, Tn m i t , 
is calculated using Equation [SI where n is the number of 
searches for a decay mode. In our search, we calculate 
90% confidence level (CL) limit, i.e., CL = 0.9 as follows: 



CL = 



r Hmi , 
r=o 



Hp(T\m)dr 



Y[P(T\ni)dr 



(8) 



r=o ; 



The lower limit of partial nucleon lifetime, t/B, is now 
calculated by: 



t/B 



1 



r 



limit 



(9) 



In the p — > vK + , prompt gamma-ray and K + — > 7T + 7r° 
search, the lower limit of the nucleon partial lifetime is 
found to be 1.0 x 10 33 years and 7.8 x 10 32 years at 90% 
CL, respectively. Only for p — > vK + , mono-energetic 
muon search, we use a special method because it has a 
lot of background. We search for an excess in the signal 
region (B3) of the muon momentum distribution. Af- 
ter applying criteria (B1,B2,B7), the events are divided 
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into three momentum bins; p^ are 200-215 MeV/c, 215- 
260 MeV/c and 260-305 MeV/c. The numbers of events 
in each momentum bin (ni, ri2, 713) are 76, 181, 185 
events, respectively. The expected numbers of neutrino 



MC events in each momentum bin (pi, 62, 63) are 78, 223, 
182 events, respectively. The nucleon decay probabilities, 
P(r|ni, n2, 713), is then calculated using EauationllOl 



P(T\ ni ,n 2 ,n 3 ) = A J J J f[- ^ + £^ + bshape > 1 * P{Y)P{\ i )P(e i )P{b)P(b shape , i )d\ i de i dbdb shape 

, 00) 

1 



where i = 1, 2, 3, corresponds to 200-215 MeV/c, 215-260 
MeV/c and 260-305 MeV/c, respectively; P{b) is defined 
as one for < b and otherwise zero; e±, £2, and £3 are 
estimated to be 0.25%, 34% and 1.3%, respectively. The 
background shape, b s hape,i, is h divided by 62- The un- 
certainty function of the background shape P{b s hape,i) 
is defined to be a Gaussian function for i = 1,3 and a 
delta function for i = 2. The uncertainties for z = 1,3 
are then estimated to be 7% and 8% from the MC model 
difference, respectively. From Equations |H1 and using 
P(r|r7 1 , n 2j n 3) instead of P(r|nj), decay limit is calcu- 
lated to be 6.4 x 10 32 years at 90% CL. Using the three 
methods, the combined lower limit on the partial lifetime 
of proton via p — > DK + is 2.3 x 10 33 years at 90% CL. 

In the n — > T>K° analysis, although we observe more 
events than the expected background in the K s — > tt + it~ 
search, there is no excess in the Kg —> 7r°7r° search. 
Therefore we also set a nucleon decay lifetime limit in 
this article. Combining two methods, we obtain the life- 
time limit of 1.3 x 10 32 years at 90% CL. In the p^ ^,+K 
and p — > e + K° search, by combining three methods, the 
lifetime lower limits are set to be 1.3 x 10 33 years and 
1.0 x 10 33 at 90% CL, respectively. The limits in all 
searches are summarized at Table [I] 

VII. CONCLUSION 

We have searched for nucleon decay via p — > DK + , 
n — > vK°, p — » and p — » e + K° from an exposure 



of 92 kt-year. No significant excess above background is 
observed. The lower limits of the partial nucleon lifetime 
at 90% CL for each mode are 2.3 x 10 33 , 1.3 x 10 32 , 
1.3 x 10 33 and 1.0 x 10 33 years, respectively (See Table[lJ|. 
From these results minimal SUSY SU(5) is fully excluded 
9J. These results also give strong constraints on other 
SUSY GUT models 0, (T3l 
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